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We study a family of frustrated anti-ferromagnetic spin-S systems with a fully dimerized ground 
state. Starting from the simplest case of the frustrated zig-zag spin ladder, we generalize the family 
to more complex geometries like tetrahedral ladders and spin tubes. After present some numerical 
results about the phase diagram of these systems, we show that the ground state is robust against 
the inclusion of weak disorder in the couplings as well as several kinds of perturbations, allowing 
to study some other interesting models as a perturbative expansion of the exact one. A discussion 
on how to determine the dimerization region in terms of quantum information estimators is also 
presented. Finally, we explore the relation of these results with a the case of the a 4 -leg spin tube 
which recently was proposed as a model for the description of the compound CU2CI4D8C4SO2, 
delimiting the region of the parameter space where this model presents dimerization in its ground 
state. 
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I. INTRODUCTION 

The Majumdar-Gosh Model (MGM)i is one of the 
paradigms in the physics of one dimensional spin chains. 
This model represents one of the first examples of sys¬ 
tems with gapped spectrum whose ground state (GS) is 
exactly known. Although the exact result is valid only 
in the point J2 = Ji/2, it results representative of an 
extended phase, a dimerized phase. A dimerized phase 
corresponds to a non-magnetic phase without “classical 
analog”, characterized by short range, strong quantum 
correlations. Systems presenting this kind of phases have 
gained interest in the quantum information community 
since, from a technological point of view, this kind of 
correlations could be exploited as a resource for quan¬ 
tum information processing^r— . On the other hand, from 
a conceptual perspective, models presenting this kind of 
non-classical phases provide a rich playground to explore 
the relations between frustration and entanglement^^—. 

Since its discovery, lots of efforts were devoted look¬ 
ing for extensions of this model to different coupling 
configuration o^dP , as well as more general lattice a^^d^ , 
and larger values of the local spin Generalizations 
of the MGM to come up with more realistic models 
are also extremely important for the theoretical descrip¬ 
tion of frustrated magnets like Cs2GuCl4i^, KGuCl^i^, 
TlGuCl^i^, NH4GuCl3i^, etc, which are being currently 
under experimental investigation. 

In this paper we present a family of anti-ferromagnetic 
spin S models with an exact dimer product state as its 
ground state. The family includes some of the aforemen¬ 
tioned models given a generalization to larger S without 
the need to include quartic or more complex many-body 
terms in the HamiltoniaiJ^. 

As in the case of the MG model, where the GS repre¬ 
sents an extended phase, the manifold in the parameter 
space where the GS can be analytically determined is rep¬ 
resentative of an extended dimerized phase, that covers 


a large region of the parameter space. 

Unlike traditional phases, which can be characterized 
in the framework of the Ginzburg Landau Theory, novel 
phases like topological or quantum spin liquids can not 
be characterized in terms of local order parameters and 
broken symmetries. In this way, quantum information 
measures have been proved to be a useful tool to charac¬ 
terize thentii, giving also information about the structure 
of the state and its correlations. For this reason, in this 
work we employ both measures of likelihood as the global 
and local fidelities to the fully dimerized state, as well as 
pairwise and block measures of quantum entanglement. 
Despite these quantities are not easily accessible from di¬ 
rect measures, they can be estimated from experimental 
parameters by measures of the structure factorsi^d^ or 
magnetic susceptibilitie s^^d^ . 

With these tools, we show that the fully dimerized 
state can be seen as the starting point to characterize 
non magnetic phases like such one observed in KGuClar^. 
This material has a zigzag structure of Cu-ions and cor¬ 
responds to a particular limit of the model studied in the 
present work. This kind of quasi-one dimensional antifer- 
romagnet structures has been intensively studie d^°i^°~ — . 

Later, we use this family as a building block to con¬ 
struct out more complex models with a similar ground 
state. In particular, we show that the ground state 
of a family of frustrated four-leg spin tubes is also a 
product state of singlets. This dimerized phase may 
appears in the strong coupling regime in three leg 
spin tubes^i. Besides, we analyze the relationship be¬ 
tween these exactly solvable models and the case of 
frustrated spin tubes, which has been recently pro¬ 
posed as a model for the magnetic behavior of the 
compound GU2GI4D8G4SO2. This compound seems 
to present frustrating anti-ferromagnetic next-nearest- 
neighbor exchange^^i^ and inelastic neutron scattering 
experiments reveal that it presents gapped and strongly 
one-dimensional excitations^!. So far, there are very few 
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theoretical studies of this frustrated model. By means of 
numerical analysis we show in this work that the ground 
state of this model presents similar features that those 
found for the exactly solvable case. 

The paper is organized as follows. In Section m we 
present the family of anti-ferromagnetic frustrated lad¬ 
ders and show that its ground state is a fully dimerized 
state. Then, a discussion about the spectrum of these 
systems and the magnetic behavior is presented. In Sec¬ 
tion m the properties of the ground state in the vicinity 
of the exact dimerizing condition are discussed in terms 
of some quantum information correlation measures. Af¬ 
terward, we show that the exactly solvable model is rep¬ 
resentative of an extended manifold. In Section IIVI the 
previous model is used as a building block for more com¬ 
plex systems which also have a fully dimerized ground 
state. For the particular case of frustrated spin tubes, 
the manifold in the phase space with partially dimerized 
GS is explored by means of numerical analysis. A nu¬ 
merical study of the ground state corresponding to the 
effective model of the compound CU2CI4D8C4SO2 and a 
comparison with the exact ground state of a very similar 
model, belonging to the family where the ground state 
can be analytically determined, are performed. Finally, 
in Section El the conclusions and some perspectives are 
presented. 


II. EXACT GROUND STATE IN SPIN-S 
LADDERS 


A. Exact manifold 


We consider the following Heisenberg model on a two 
legs spin-S' ladder 

N 

H = ^ J{i) S 2.-1 • S 2 . + J'{i) S 2 . • S2(,+i)_i 

+ J" {i) S 2 i_l • S 2 (i-|- 1 ) -I- J 2 {i) §21 • S 2 (i+ 1 ) ( 1 ) 

+ <^2(0 S 2 i -1 ■ S 2 (i+i)_i , 

where represents the local spin on the site k and = 
Sfc_|_2Ar, with N being the number of rungs. This ladder is 
represented in Figure[T]-A. Let us consider the case where 
all the couplings are positive (all the interactions are anti¬ 
ferromagnetic). Starting from the general spin-S model 
without translational invariance in Eq. 0 we can show 
that, imposing a simple constraint on the couplings on 
each square plaquette defined by the sites { 2 i—l, 2 i, 2 (i-|- 
1 ) —l, 2 (i-|-l)}, the ground state of the system is the fully 
dimerized state lijj) = {^^4 |0)i, with 


|0)z = 


y/ 2 S+l 


E(-i) 


m+51 


( 2 ) 


m— — S 


where index i labels the rung in the ladder and 
|m, —m)i are product states such that S2i_i\m, —m)i = 


— S2j|m, —m)i = m|m, —m)i on the rung i. In order 
to show that IV') results an eigenstate of H with energy 
Eq = —JNS{S + 1 ), we rewrite the Hamiltonian in terms 
of local operators on each rung 


Li — S2i + S2i_l (3) 

K* = S 2 . - S2,_i. (4) 


Here, is the total angular momentum of the rung i, 
and Ki is a set of local observables which completes the 
full local Lie algebra of observables: 


[L^i, Ljy] 


( 5 a) 



( 5 b) 

:k^,k,] 


( 5 c) 


where is the fully antisymmetric Levi-Civita symbol 
and i is the imaginary unit (i^ = — 1 ). In terms of these 
rung operators the Hamiltonian reads 


N 


i=l 

N 




+ J"(z) 


N 

+ E 

2=1 


(k. 


• Li+i — Li • Ki-i-i 


If J'{i) + J''{i) = 2 J 2 {i), the last term in ([ 61 ) vanishes and 
the state \tjj) is an eigenstate of the Hamiltonian due to 
Li|0)i = 0. Noteworthy, this result is valid for any value 
of the local spin magnitude S. 

Now, provided the condition J{i) > (z — 1 ) -|- 

J {i) + J (z — 1 ) -I- J"(z)), |z/’), we show that it is the only 
ground state of the system. For this purpose, we rewrite 
the Hamiltonian in a convenient form (see Figure [ 1 ]): 
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FIG. 1 . (Color On-line) Tetrahedral-ladder geometry and 
its decomposition in terms of three Hamiltonians with fully 
dimerized GS. 


H = Hs -f He + Ho + 

1 ^ r - 

Ho = - J (f — 1) ^S2i-2 + S2i-1 + S2i^ 

^ i=l L 

-\-J (i) ^S2i-l + S2i + S2i-|-1^ + 

+ {j'{i - 1) + j'{t))S {S2^-l + S2,f- 


He = 


N 


ic = ^ ^ /'(* - 1) (S2.-3 + S2.-1 + S 

2 


2 i 


+ J (i) ^S2i-l + S2i -|- S2i-|-2^ 

+ (/'(! - 1) + J"(z))^ (S2._i + S 2 O'- 
-(j"(f-l) +j"(z))B(B+l) 


N 


Ho = E 


A(*) 


(S2.-1+S2.)' 


(7) 


( 8 ) 


(9) 

( 10 ) 


where Ed = —S{S -|- 1) J{i) is the energy for the 


state and 

A(i) = J{i) - (i - 1) + /(*) + /' (*-!) + 

is an effective coupling constant associated to L^, the to¬ 
tal angular momentum of the pair. Ho^e correspond, up 
to a constant, to the Hamiltonians of two zig-zag ladders 
(see Fig. [T]-B and[T]-C), while Ho is the Hamiltonian of a 
set of uncoupled pairs (Fig. [IJ-D). Now we will show that 
Ho and He are semi-definite positive operators, being 
|'0) its ground state. 

To see it, we observe that the minimum eigenvalue of 
a sum of operators is bounded from bellow by the sum 
of the minimum eigenvalues of its terms: 

min A > 7 min Ai, 

AeA(5:ihi) ^AieA(hi) 

where A( 0 ) = {Ai, A2,...} denotes the spectrum of op¬ 
erator O. Now, we notice that for J'{i — l),J'{i) > 0 , 
due to the theorem of addition of angular momentum, 
each term in ([ 5 ]) and (0 is bounded from bellow by 

Ai > b ~ b + J b) ^s-ii\i\s-ii\+i)-sis+i)-siii+i)ii = 0 

4 = 

and hence, Hb,Hc > 0 . Here, k is the total spin in 
rung i. 

On the other hand, is easy to verify that Hb|' 0 ) = 0 , 
from which it follows that |'0) is an eigenvector of Hs^c" 
with minimum eigenvalue. Since for each i, J(i) > 
— + J {i) + J (z —1)-I-J"(i)), Hu also results 

positive, and hence, H > Ed- But \'ip) is an eigenstate 
of H which saturates that bound, so it is a ground state 
of H. Due Hu is gapped, \ijj) is the unique state that 
saturates the bound. 

For the S = 1/2 case, we can improve this bound 
by observing that minA((S2i-i -I- 82^ -|- S2i-i-i)^) = 
minA((S2i-2 + 82*-! + 82i)^) = 3 / 4 , disregarding the 
value of (82^-1 -t- 82^)^. This allows us to move the terms 
in (82i-i+82i)^ in Hu and Hq to Hu, leading to the im¬ 
proved bound J(z) > |(J (z— 1 )-|-J {i)+J {i—l)+J''{i)). 

Finally, for the translational invariant case, J{i) = J, 
J'{i) = J', J"{i) = J", the sufficient condition for the 
exact dimerization is given by J' + J" = 2 J2 and 

J(z) > I ( 11 ) 

^ ^ (S+1)(U + J//) S>i/2 

The Majumdar-Ghosh point can be recovered from this 
result as a limit. For S = ^, homogeneous couplings 
J ^ J' = 2J2 and J” = 0 the GS is still a dimerized 
state but it is degenerate. 


B. Localized triplons 

An interesting subfamily of models corresponds to the 
symmetric Tetrahedral ladder, which is obtained by set¬ 
ting J'{i) = J"{i) = In this subfamily, 
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FIG. 2. (Color On-line) Ground State Energy per bond cor¬ 
responding to the S = 1/2 ladder with J' = J" (blue on-line) 
and the Ground State Energy corresponding to a S' = 1 chain 
(red on-line) as a function of J 2 calculated with DMRG for 
N=60 rungs(sites). In the blue light region the ladder has 
a dimerized ground state whereas in the orange region the 
S = 1/2 ladder is equivalent to the S = 1 Heisenberg chain. 



those terms in m containing operators K vanish and 
the Hamiltonian depends only on the total spin of each 
rung Li. 

N 

H-^J(*) 

2 = 1 

Since [Lf,Lj] = 0, [H, Lf] = 0. Therefore, each 
eigenspace of H can be characterized by {k}, the 
set of total angular momentum quantum numbers 
(Lf |{Zi},...) = l{l + l)|{/i},...)) associated to each 
rung. Hence, on each proper subspace, the model re¬ 
duces to a spin chain with different values of the spin at 
each site {k = 0,1,2...2S) with on-site quadratic terms 

and exchange terms Li ■ L^+i. If the total 
spin in a given rung is zero, there is no coupling with its 
neighboring rungs. 

For S' = 1/2 the ground state of the system may 
correspond to = 0 or = 1, this is, the elemen¬ 
tary excitations of the system can be seen as localized 
triplons. The number of these triplons in the ground 
state is determined by the competition between the terms 

J{i)^ and J 2 {i)Li ■ L^+i. By setting the condition 
(S -I- 1)(J'(*) -I- J'{i — 1)) < J{i), IV') is the ground state, 
corresponding to {k = 0}. On the other hand, large val¬ 
ues of J 2 favors larger values of k. In Figure [2] we show 
the energy per bond for a S = 1/2 ladder as a function 
of J 2 with J' = J" (blue circles), calculated by means 
DMRG^. Red circles corresponds to the energy of the 
GS in the sector li = \. It is clear that for J 2 > 0.7 the 
spin ladder behaves like a S' = 1 spin chai n^^d^iSO , 

This result is important for two reasons. On the one 
hand, it gives us a picture about how the dimerization 
breaks when we cross the boundaries of the exact mani¬ 
fold: for large enough J', the system suffers a level cross¬ 
ing to a higher value of the rung local spin k. On the 


(^^-S(S+1)^ +J'(z)L, •£,+!. 


( 12 ) 


other hand, we can take advantage of this result to gain 
information about the magnetic behavior of the system. 

We introduce a coupling term to an uniform external 
magnetic field h: 

+H. (13) 

i 

For the S = 1/2 case the magnetic process has been 
studied by mean a strong coupling approach^. The ex¬ 
ternal magnetic field induces a competition between the 
one-site term and the magnetic field term. If all J{i) are 
large, the GS remains being |V’) up to a magnetic field 
value (lets say hi) where there is a crossing level and the 
GS becomes the state with singlets in the even/odd rungs 
and fully polarized spin-1 states in the odd/even rungs. 
Notice that for these states the exchange interaction does 
not contributes. 

Then, if we increase more the magnetic field, we ob¬ 
tain a second crossing level with the state containing 
fully polarized spin-1 states in all the rungs. In the 
same way, we will find successive crossing levels with 
fully polarized states with total spin lodd and leven in 
odd and even rungs. Each one of these crossing levels 
result in a jump in the magnetization curve followed by 
a magnetization plateau. Hence, the magnetization pro¬ 
cess is given by successive jumps and plateaux. When 
we slightly move away from the condition J'{i) = J"(i), 
these plateaux may remain present in the magnetization 
curve, whereas the jumps are smoothed. This frustration 
induced plateaux has been studied forS'=l/2,S'=l 
and S' = 3/2 ladders^^— The present analysis provides 
a simple theoretical explanation for this behavior. 


III. VICINITY OF THE EXACT MANIFOLD 

Although we have shown that the state |V’) is the GS 
of the model ([6]) just when the condition J'{i) -|- J"(i) = 
2 J 2 (I) is fulfilled, this state is representative of a region in 
the parameter space, where results an accurate approxi¬ 
mation to the true ground state. In order to character¬ 
ize the region presenting dimer order, we consider both 
measures of similarity between Itjj) and the GS, as well 
as measures of entanglement. 

The Uhlmann’s Quantum Fidelity^ provides a measure 
of the similarity between two quantum states: 


= Try ^/ap^/a = F[<j, p] 

where p and a are two quantum states of the same sys¬ 
tem. If a corresponds to a pure state |a), this quantity 
reduces to F[p, |a)] = {a\p\Q). 

We start considering the fidelity between |^) and the 
true GS in a region close enough to the exact manifold 
such that a first order perturbative treatment would be 
feasible. Starting from the Hamiltonian ([HI), and using 
the algebraic properties of and the canonical 
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first order perturbation theory leads to 


|GS)«.(1-A7V^2)1/2|^)^1^ 


7*lb* + l), (14) 


where \i,j) = • K^-|' 0 ), a 

'^’('S’+i) J'h)+J"(b- 2 J 2 (i) J 2 _ J_ ^2 Notice 

the explicit SU (2) invariance of the approximation, 
as well as its translational invariance for the homoge¬ 
neous case. In this way, the global fidelity is given by 
(1 — which is valid for <C 1. From 

now on, we are going to restrict to this last case and 
hence, = ± 7 . In this case, the previous result seems 
suggest that for large systems the dimerization is con¬ 
strained just over the exact manifold. However, to look 
for high values of the fidelity in a large system is a very 
demanding condition. Due to its definition, for product 
states ^"[ 10 !)®^, |a')®^] = (J^[|a), |a')j)^ and hence, yet 
for very similar states, the fidelity vanishes in the large 
N limit. On the other hand, from (fT^ we can estimate 
the fidelity for the state of a single rung (P 12 ) against the 
singlet state J'q[pi 2 ] = J^[pi 2 , |singlet)] Ri (1 - 
When this approximation is valid, we can see that 
J^[|GS), 1^/;)] ~ To[pi 2 ]^, which is just what we expect 
if the global state behaves like a product of the local 
states of the rungs, which is an important feature of 
the dimerized phase. A similar result can be obtained 
by the method of variational cluster mean field -I- RPA 
discussed in^. This treatment predicts that the GS is 
well approximate by |^/)) plus small Gaussian correlations 
for I 7 I < 1 . 

Now, we will see that for the S' = 1/2 case, the value 
of J^o[pi 2 ] determines most of the relevant features of the 
dimerized phase. To see this, we observe that due to the 
SU{2) symmetry, for the S = 1/2 case the local state pij 
of a subsystem composite by the (single spin) sites (z,j) 
is completely determined by Fo[pij\-. 


1 

Pij = - 5 - 14 


4-Ao[Pij] 1 


(15) 


Since pij is the state of a subsystem associated to pure 
global state, a measure of the correlations between this 
subsystem and the rest of the system is given by its en¬ 
tanglement entropy^. For the state dTSl) it is reduced to 


S{pij) = h {J^oiPij]) + 3/1 


1 -^olPij] 


where h{x) = —a;log 2 (a;). For the dimerized phase, 
where the global state is well approximated by the fully 
dimerized state, which is a product state, S{pi 2 ) should 
remain small, while S{pij) ~ 2 for j ^ 2. In order to give 
a more accurate idea about the limit value of -Fo[pi 2 ] for 
which it makes sense to talk about dimerization, it is 
better to analyze the internal degree of entanglement of 
the rung. A measure of the degree of pairwise entan¬ 
glement for mixed states is provided by the logarithmic 



FIG. 3. (Color On-line) Relation between entanglement and 
fidelity to the singlet state for SU (2) invariant states of pairs 
of 5" = 1/2 systems. For J-o[p\ < ll\/2 the internal degrees of 
freedom are not entangled. 


negativity^^*^ 

^A\b[p] = log2 

where |A|i = tiV .A is the trace norm (or the Schatten 
— 1 norm) of the matrix A, and represents the par¬ 
tial transposition of p with respect to the subsystem A, 
i.e. the linear map defined by (|Q:)^|/3)B(a'|^(/I'|B)*'^ = 
|a')^|/3)B(Q!U(/3'|B- This quantity is saturated by 
f-™ [jsinglet)] = log 2 (l -I- 25') and vanishes for every sep¬ 
arable state. For the state is given by 

£^j[Pij] = log2(l + max(0,2J'^[pij] - 1)). 

In Figure [3l the behavior of 5[pi2] and E^^[pi 2 ] as a 
function of J^o[Pi 2 ] is depicted. Notice that for Jo[Pi 2 ] < 
the reduced state is separable and hence, the corre¬ 
spondent global state is not dimerized anymore. 

In Figure m a landscape of the fidelity between the sin¬ 
glet state and the state of the rung J'o[pi 2 ] (A) and a 
nearest neighbour external pair J^o[p 23 ] (B), obtained by 
numerical evaluation is shown, for the case of the sym¬ 
metric ladder (J" = J'). The presented results were eval¬ 
uated by means of the Lanczos method^^. The straight 
dashed line (red online) indicates the intersection with 
the exactly dimerized manifold where the fully dimerized 
state is the ground state of the system. For the strong 
pair (A), the fidelity is symmetric regarding the exchange 
between J' and J 2 . The graphic reveals a wide dimerized 
region (Fo[pi 2 ] > 0.95) around the exact line J' = J 2 , im¬ 
plying that the pertubative analysis is accurate over this 
region. As we cross the critical value J' ~ 0.6J, Fb[pi 2 ] 
is suddenly reduced, due to the GS is now orthogonal 
to the dimerized one. For pairs of spins coupled by J' 
(B), we observe that near the exact line Fb[P 23 ] ~ 1/2, 
which is consistent with a fully mixed state. For larger 
J' these pairs become more entangled, at the expense of 













6 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 


FIG. 4. (Color On-line) Landscape of the fidelity between 
the singlet state and the state of the internal pair pi 2 (A) and 
state of the nearest neighbor external pair p 23 (B) associated 
to the gronnd state for the spin 1/2 ladder with J' = J" for 
N = 8 rnngs. The straight dashed line (red on-line) indicates 
the exact dimerized line J 2 = J^ In the panel (A), the dark 
continnous line indicates the boundary of the region with fi¬ 
delity > 0.95. In both panels, the thick dashed cnrve (yellow 
on-line) corresponds to the valne l/%/2, limiting the region 
where the correspondent local states are entangled. Insets 
indicate the pair used to calculate the fidelity. 


the entanglement in the rung (1 — 2). On the other hand, 
increasing J 2 the spins on the pair 1 — 3 tends to align, 
which reduces its fidelity to the singlet below 1/2. 

Further improvements can be obtained by mean of a 
cluster mean field -I- RPA expansion as we have shown in 
a previous worb^i. This opens a way to solve some effec¬ 
tive models associated to certain quasi-one dimensional 
materials like KCuCla. 



FIG. 5. (Color On-line) Building a spin tube family with a 
fully dimerized ground state. 

IV. FRUSTRATED FOUR-LEG SPIN TUBES. 

The family of ladders presented above can be used as a 
“building block” to obtain more complex models with a 
product state as its ground state. It is straightforward to 
show that there is a family of Hamiltonians correspond¬ 
ing to four-leg spin tubes that present also a dimerized 
ground state. These Hamiltonians can be written as a 
sum of two ladders whose Hamiltonians are given by 0 
and N square plaquettes as represented in Figure 0 

N 

^^tube — H/fjdderl ^iladder2 ^ ^ H square^j (16) 

i=i 

Where N is the number of rungs in the ladders. Hiadderi 
and Uiadder 2 are Hamiltonians corresponding to the up¬ 
per and lower ladders in Fig. [5] and Usquarej is the 
Hamiltonian of the j-th square plaquette in the fig¬ 
ure. These square plaquettes have also a dimerized 
ground state corresponding to dimers in the strongest 
bonds, as can be easily seen since the square is a spe¬ 
cial case of ladders discussed previously with 2 rungs. 
Let us concentrate in the homogeneous case J(i) = J, 
J'(z) = J' and = J", Vj, keeping in mind that 

all the conclusions can be easily generalized for the in- 
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homogeneous case. Different combinations of couplings 
J', J" and J 2 preserving the constraint J' + J" = 2 J 2 
gives different geometries of four-leg spin tubes with 
a fully dimerized ground state that can be written as 
\tptube) = |V');adderix|'0)iadder'2- Then, following the Same 
steps as in the previous section is easy to show that the 
state I'(/’tube) is an eigenstate of the system and there is 
a range of couplings where this state is the ground state 
of the system. Besides, as in the ladder case, we can ex¬ 
pect this exact GS to be representative of a finite region 
around the exact manifold of the parameter space. 

We will consider now the particular cases of tubes with 
J' = J” and J" = 0. For these two cases, the fully 
dimerized state \il)tube) is the ground state of the system 
if the following condition is satisfied 


J' < 


S=\ 

^(s+i) S >\ 


(17) 


where ^ = 1/3 for J' = J" and ^ = 2/3 for J" = 0. 
In Figure |6] the energy per rung in units of S{S -I- 1) as 
a function of the coupling J 2 , corresponding to the case 
J" = 0 and J' = 2 J 2 is depected for different values of 
the spin. The range of values where the ground state is 
the product singlet state {i.e. E/N^ungs = —1) is bigger 
than that found analytically. The reason is that Eq. uni) 
represents just a sufficient condition, being the real range 
larger in general. 

As we have seen for the ladders, the symmetrical case 
J' = J" is special. The excited states correspond to lo¬ 
calized triplons and the eigenvalues of the Hamiltonian 
can be labeled by the set of values of the total momentum 
in the rungs {h}- This situation gives a rich magnetiza¬ 
tion profile, containing a sequence of jumps and plateaux. 
Such scenario is due by the competition between the 
terms in the Hamiltonian that favors states with minimal 
total momentum in each rung, spin-exchange terms and 
the magnetic field contribution. The crossover between 
the ground state \iptube) and the state with alternation of 
one singlet and one triplet in each rung determines the 
jump in the magnetization curve between the plateaux 
at m = 0 and m = 

In this case, system frustration promotes the magne¬ 
tization plateaux and the jumps in the magnetization 
curve but also makes the ground state simpler (in each 
magnetization sector the ground state is a direct product 
of rung-states). 

Even while the condition J' + J" = 2 J 2 is hold, for 
J" J' this is not true any more. Changing the value of 
J"/ J' the magnetization plateaux reduce its widths and 
jumps between plateaux transform into a smooth piece 
of the magnetization curve. 

The existence of the exact result for tubes also leads 
us to ask in which conditions, two ladders which present 
dimer order in its ground state, conserve it when they 
become weakly coupled, assembling a tube. From a per¬ 
turbative argument, for small enough inter-ladder cou¬ 
plings, we expect that the global state stays dimerized. 



FIG. 6. (Color On-line) Ground state energy per bond in 
units of S{S -I- 1) as a function of J 2 calculated with DMRG 
for a four-leg spin tube with 160 spins, J” = 0 and J' = 2 J 2 . 
Blue circles, red squares and yellow rhomb! correspond to 
S = 1/2, 5 = 1 and S = 3/2 respectively. 


However, due to the exact result Eq. ca), despite it be¬ 
comes more frustrated, for larger couplings the system 
approach to another exact dimerized configuration. As 
a result, the region presenting dimer order would be en¬ 
larged. As an example, we will consider the case in which 
the interaction between ladders are given between corre¬ 
spondent spins on each ladder and on one of the diagonals 
of each plaquette (see the inset of Figure [7]), calling Jy 
and Jd the respective coupling constants. This case is in¬ 
teresting since recently, a similar topology was proposed 
(but for Jd = 0) as the appropriate model describing the 
compound Cu2Cl4D8C4S02^i^. Despite this case does 
not satisfy the exact dimerization condition for tubes, it 
is interesting to find out if it could support dimer order. 
In Figure [7] the behavior of the fidelity between the state 
of a rung and the singlet state, as well the entanglement 
entropy of this subsystem with the rest of the tube is de¬ 
picted, for weakly coupled zig-zag ladders, as a function 
of the inter-chain couplings. Notice that the dimeriza¬ 
tion over the lateral ladders is not broken for quite large 
values of Jy and Jd near the exact dimerization condi¬ 
tion. It would suggest that we can expect the presence of 
dimer order in the model proposed for CU 2 CI 4 D 8 C 4 SO 2 . 
In the next sections we extend this result for a more re¬ 
alistic case, when the ladders does not satisfies the exact 
dimerizing condition. For it, we take advantage of the 
exact result for the four legs frustrated tube discussed 
above starting from a highly frustrated system but with 
a separable GS. 


A. Effective Hamiltonian of the four-leg spin tube 
material Cu2Cl4DgC4S02 


In a Recent experiment, inelastic neutron scattering 
has been used to investigate the magnetic excitations in 
the quantum spin-liquid system Gu2Cl4D8G4S0 2^^'^^ . In 
that work, it was suggested that the appropriate Heisen- 










Jd 


FIG. 7. (Color On-line) Dimerization in the case of weakly 
coupled ladders as a function of inter-ladder conplings calcu¬ 
lated by using the Lanczos method with L = 4, J 2 = 0.5J' 
and J' — 0.95 J (see the inset) . (A) Entanglement entropy of 
a strong coupled pair with the rest of the tube. (B) Fidelity of 
the local state with the singlet state. The region in which the 
fidelity is over 0.99 matches with such that the entanglement 
entropy is lower that 0.2. The dashed thick line (red on-line) 
corresponds to the exact dimerization condition J± = 2Jd 



els are similar, the effective model for CU 2 CI 4 D 8 C 4 SO 2 
may be studied starting from the exactly known ground 
state, taking the diagonal couplings in the squares as a 
perturbation. Although this perturbative study is out 
of the scope of the present paper we can see that the 
dimerized state is robust in the exact model. This ro¬ 
bustness suggest to make an expansion around the dimer 
ground and incorporating triplon excitations. In the rest 
of this section, we are going to analyze numerically the 
proposed model for the material CU 2 CI 4 D 8 C 4 SO 2 , look¬ 
ing for fingerprints on the properties predicted for the 
exactly solvable case, leaving the analytical study of the 
corrections for a future work. 


B. The model for the material and the exactly 
solvable model 



berg Hamiltonian is a S' = 1/2 four-leg spin-tube with 
no bond alternation as the showed in Figure [ 8 ]-a). There 
are a scarce number of theoretical results on this kind 
of prototypical models on the spin tube a^^i^^ . We study 
numerically the frustrated four leg spin tube model pro¬ 
posed to describe the compound CU 2 CI 4 D 8 C 4 SO 2 . 

The proposed model is closely related to the family 
of spin tubes presented in the previous sections. Con¬ 
sider the member of the family of spin tubes with a fully 
dimerized ground state schematized in Figure [ 8 ]-b). This 
model can be obtained from the Hamiltonian proposed 
for the material adding an extra diagonal coupling in 
each square. This modified Hamiltonian belongs to the 
family of spin tubes presenting a dimer product ground 
state. If the ground state properties of these two mod- 


Above we have analyzed the behavior of the ground 
state in tubes near the dimerization condition. Now, we 
are going to discus which are the common features be¬ 
tween these results and the behavior of a more realistic 
family of tubes. In particular, we consider the model pro¬ 
posed by Garlea et al. i n^^'^^ for the CU 2 CI 4 D 8 C 4 SO 2 
compound. In order to see what features are shared be¬ 
tween our exactly solvable model and the model proposed 
i n25,26 ^ we explore numerically the ground state of Hamil¬ 
tonians on the family J > J 2 , J', corresponding to the 
tube in the figure [5]-a, looking for common features. 

In Figure [9] a landscape of several entanglement ob¬ 
servables, corresponding to the ground state of a tube 
with couplings as in Figure [5]-a are depicted. In the top 
panels, the fidelity between the local state of a single 








































9 



FIG. 9. (Color On-line) Local structnre of the ground state for the tube, as a function of couplings J 2 and J', evaluated through 
exact diagonalization (Lanczos method) for the case L = 4 plaquetes (16 sites) and S = 1/2. Panel A: Fidelity of the local 
state associated to a single rung respect to the singlet state. Panel B: entanglement entropy of a rung with the rest of the 
tube. Panel C: Entanglement of a plaquette; Panel D: internal entanglement, measured by its logarithmic negativity between 
two pairs in the same plaquette. The dashed (red on-line) straight line represents the condition J' = 2 J 2 ; the dark continuous 
curve bounds the region where J^o[pi 2 ] > 0.95 while the dashed light curve (yellow on-line) is the bounds the region where the 
pair pi 2 is entangled (J^o[pi 2 ] > l/%/2). 


rung and the singlet state (Panel A) and the entangle¬ 
ment entropy of a rung (Panel B) are shown. Due to 
the SU (2) symmetry, for local states of a single rung the 
entanglement entropy is a function of the fidelity to the 
singlet state. As we could expect from a composite mean 


field treatment^i^, the local state of a single rung can 
be accurately approximated by a pure singlet state in 
a relatively wide region around the condition J' = 2 J 2 , 
even if J 2 is moderately large. 

On the bottom panels we can appreciate the entangle- 
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merit of a plaquette composed by two parallel rungs with 
the rest of the tube (Panel C) and the internal entangle¬ 
ment between two rungs in the same plaquette, measure 
by its logarithmic negativityS^ (Panel D). For large J 2 and 
J', we observe that both the entanglement between the 
state of a square and the rest of the tube becomes larger, 
which is compatible with a symmetry broken phase. In 
fact, the limit J —>■ 0 with J' and J 2 fixed corresponds 
to a non frustrated tube, for which a Neel-like phase is 
expected. We can also notice that the region where the 
squares are not entangled is larger than the dimerized 
region. In particular, near J 2 ~ 0.5J and J' ~ 0.9J, 
we can observe a region where the entanglement of the 
square with the rest of the tube is small, but the inter¬ 
nal entanglement between rungs is near to 1. This seems 
to indicate that such region corresponds to a resonant 
plaquette order, where the GS is well approximated by 
\GS) Ri |a)^„Qre, being \a)square a four spin singlet state 
i^squarel^) square = 0 )j which is not a product state of the 
rung states \a)square |single)]singlet). 

In summary, we saw that although the exact dimer¬ 
ization condition for this model is only possible in the 
trivial limit J' = J 2 = 0 , the phase diagram looks quite 
similar to that associated to the family containing the 
solvable model. This can be understood by considering 
the removal of the extra bond on the square hamiltonian 
as a perturbation over the Hamiltonian of the solvable 
tube. However, near the crossover point, the perturba¬ 
tion theory is no longer valid, which can give place to 
new features in the phase diagram. 

V. DISCUSSION AND PERSPECTIVES 

In the present paper, a general SU (2) invariant quan¬ 
tum spin-S Heisenberg ladder was investigated. A suffi¬ 
cient condition for the existence of a fully dimerized exact 
ground state was shown for a wide subfamily of such sys¬ 
tems. Besides, by means of a combination of numerical 
and analytical techniques, the existence of this phase for 
a general value of the local spin was proved, showing that 
the region in the parameter space corresponding to the 
dimerized phase is reduced as the magnitude of the local 
spin grows. 

For the case of symmetrical frustration the excitations 
were also exactly determined and a discussion about the 
magnetization process^Sr— was presented. 

The ground state properties around the exact manifold 


in the parameter space were explored by means of numer¬ 
ical analysis, showing that these remain close to the exact 
case over a finite region. Due to the large stability of this 
phase against external perturbations, a quantum simula¬ 
tor that could reproduce this kind of couplings would be 
able to prepare a large number of fully entangled pairs 
in a robust way. In the last years several proposals for 
the experimental simulations of spin systems in ion trap 
experiments^^— suggest that this kind of setup could be 
readily in a near future. 

Besides, we have shown how the family of Hamiltoni¬ 
ans with a fully dimerized ground state can be extended 
from the family of ladders to more complex models. As 
an example a family of frustrated 4-leg spin tubes with 
a dimerized ground state was built. Afterward, common 
features found in the ground state of the solvable family 
and those obtained for more realistic models were ana¬ 
lyzed. In particular, a comparison to the model proposed 
for the CU 2 CI 4 D 8 C 4 SO 2 compoun d^^i^^ was discussed. 

We hope this study can be taken as a starting point for 
a more systematic study of the ground state properties 
and excitations around the lines on the parameter space 
where the ground state was exactly determined. 

As a perspective, the study of hole dynamics on a back¬ 
ground of dimers in one-dimensional systemsii. In higher 
dimensions there is an important amount of results on 
dimerized ground states^^— that can be used as start¬ 
ing point to study hole dynamics. In 2D has been re¬ 
cently proved that the quantum statistics of holes in a 
dimer background can be changed without affect the en¬ 
ergy dispersion^ and the density of holes has an impact 
on the magnetization plateaux^. As we start from these 
families of models where the ground state is a dimer cov¬ 
ering, introducing holes in the system may result in a 
very interesting phase diagram. 
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